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Lorentz symmetry is one of the cornerstones of mod¬ 
ern physics. However, a number of theories aiming at 
unifying gravity with the other fundamental interac¬ 
tions including string field theory suggest violation of 
Lorentz symmetry (THU. While the energy scale of such 
strongly Lorentz symmetry-violating physics is much 
higher than that currently attainable by particle ac¬ 
celerators, Lorentz violation may nevertheless be de¬ 
tectable via precision measurements at low energies El. 
Here, we carry out a systematic theoretical investiga¬ 
tion of the sensitivity of a wide range of atomic systems 
to violation of local Lorentz invariance (LLI). Aim of 
these studies is to identify which atom shows the biggest 
promise to detect violation of Lorentz symmetry. We 
identify the Yb + ion as an ideal system with high sen¬ 
sitivity as well as excellent experimental controllability. 
By applying quantum information inspired technology 
to Yb + , we expect tests of LLI violating physics in the 
electron-photon sector to reach levels of 10 -23 , five or¬ 
ders of magnitude more sensitive than the current best 
bounds (5H3. Most importantly, the projected sensitiv¬ 
ity of lO -23 for the Yb + ion tests will allow for the first 
time to probe whether Lorentz violation is minimally 
suppressed at low energies for photons and electrons. 

Formally, we can classify LLI-violating effects in the 
framework of the Standard Model Extension (sme) 03. 
The SME is an effective field theory that maintains Lorentz 
invariance of the total action, energy-momentum conserva¬ 
tion, and gauge invariance, but supplements the Standard 
Model Lagrangian with all combinations of the SM fields 
that are not term-by-term Lorentz invariant. Here we focus 
on the c^y tensor term of the SME Lagrangian signifying 
the dependency of the maximally attainable velocity of a 
particle with respect to its propagation direction. SME al¬ 
lows for a violation of LLI for each type of particle, mak¬ 
ing it is essential to verify LLI in different systems at a high 
level of precision. As a result, LLI tests have been con¬ 
ducted for the photons 0 , protons 0, neutrons (9j 10], 
electrons 00, and neutrinos m with the detailed sum¬ 
mary of all current limits given in ff2l . 

Testing LLI of the electron motion in a Coulomb po¬ 
tential created by a nucleus has the appeal of testing for 
new physics in a well understood system. In these atomic 
experiments 00 , one searches for the variations of the 
atomic energy levels when the orientation of the electronic 
wave function is rotated with respect to the hypothetical 
preferred reference frame. The analysis of any experiment 


(a) Lorentz violation: 
electron sector 



(b) Lorentz violation: 
photon sector 



FIG. 1: To analyze the atomic LLI experiments, we need 
to pick a reference frame which allows for two 
interpretations of the result: (a) we can either assume that 
the Coulomb potential is symmetric and any 
Lorentz-violating (LV) signal is attributed to the electron; 
or (b) we assume that electron obeys Lorentz symmetry 
and any LV signal is attributed to the photon sector. 


requires the selection of the preferred reference frame. In 
the present context, this lends itself to the interpretation of 
the results in terms of the Lorentz symmetry test in either 
electron or photon sectors as illustrated in Fig. |T] 

In the present work, we discuss how to drastically 
improve LLI tests in the electron-photon sector. We 
identify several factors that affect the precision of the LLI 
test with atomic systems: 

1. The lifetime of an atomic state used to probe Lorentz- 
violating effects. 

2. The sensitivity of this state to the LLI-violation effects, 
i.e., the size of the matrix element of the corresponding 
operator. The LLI violating interaction 0 contains 
second power of the momentum operator p. Therefore, 
corresponding matrix elements are expected to be large 
for orbitals with large kinetic energy. This happens for the 
atomic 4/-electrons localized deep inside the atom in the 
area of large potential and kinetic energy in some atomic 
systems. 

3. Good theoretical understanding of the system and 
availability of already developed experimental techniques. 

We find that the metastable 4f u 6s 2 2 F 7/2 state of Yb + 
has extremely favorable characteristics such as being very 
sensitive to LLI violations as well has having an exception- 
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ally long lifetime on the order of 10 years m. In addi¬ 
tion, the experimental methods for Yb + are well developed 
in the context of atomic clocks and measuring the varia¬ 
tion of fundamental constants. In particular, the electric- 
octupole E3 transition between the 4/ 13 6s 2 2 F 7 / 2 excited 
state and the ground state 116 -181 has been studied in de¬ 
tail for these applications. Yb + ions are also used in quan¬ 
tum information research oa. As a result, the Yb + ion is 
a well-studied system for precision control and manipula¬ 
tion of its atomic states making it particularly well-suited 
for searches of Lorentz violation signature. 

The relevant electronic Lagrangian in the SME (in 
atomic units) is given by 

z = l;<A(y v + c'^y) D v ip-pip, (1) 

where if/ is a Dirac spinor, y^ are the Dirac matrices, 

if/ D v if/ = if/D v if/ - il/D v if/ with D v being the covariant 
derivative. The tensor + k^/2 that character¬ 

izes the LLI-violation effects contains Lorentz-violation 
parameters from both the electron (c^y) and the photon 
(k^y) sectors mm as illustrated by Fig. [I] 

From Eq. Q. violations of Lorentz invariance and Ein¬ 
stein’s equivalence principle in bound electronic states re¬ 
sult in a small shift of the energy levels described by a 
Hamiltonian l6l 

where p is the momentum of a bound electron, c is the 
speed of light, and U is the Newtonian gravitational po¬ 
tential. The parameters cf\ coo, and C (2) are elements 
in the c ^ tensor. The relativistic form of the p 2 opera¬ 
tor is cyoyipj. The nonrelativistic form of the 7^ 2) op¬ 
erator is Tq 2) = p 2 - 3p 2 , and the relativistic form is 
Tf = cyo (y j Pj - 3y 3 p^, with z (the 3rd spatial compo¬ 
nent) defined by the quantization axis. 

The value of c^y is specified in the Sun’s rest frame. Be¬ 
cause of the Earth’s motion, time dependent-Lorentz trans¬ 
formations from the Sun’s rest frame to the local labora¬ 
tory frame give rise to time-dependence of the local ob¬ 
servables that involve any of the C[ } 0) , coo, and/or C^ 2) pa¬ 
rameters. The elements cjk in c^ y where /, K = {1,2,3} 
(spatial components), which describe the dependence of 
the kinetic energy on the direction of the momentum, have 
a leading order time-modulation period related to the side¬ 
real day (12-hr and 24-hr modulation). Other elements, 
ctj (where T = 0) and coo, describe the dependence of the 
kinetic energy on the boost of the laboratory frame, and 
have the leading order time-modulation period related to 
the sidereal year. 

Yb + ground state has [Xe]4/ 14 6s ground state, and is 
essentially a heavier analog of Ca + which was used in the 
most recent experiment □ , but has additional set of states 
resulting from excitation of the filled 4/ shell. Yb + also 
has a metastable excited 5ds/2 state analogue to Ca + . How¬ 
ever, the 4/ shell in Yb + is localized much deeper in the 
atom than the 6s or any other valence states, such as 5 d. 



^ Cr 



FIG. 2: In a magnetic field B, two orthogonal orientations 
of the electronic wave function of the 4/ 13 6s 2 2 F-U 2 mj 
manifold in Yb + , with mj = |1/2| and mj - |7/2|, will 
have a different response to the Lorentz-violating effect 
quantified by the c^ tensor. This results in the energy 
difference between these two states while the Earth rotates 
the wave function with respect to the hypothetical fixed 
reference frame. 


As a result, we expect the 4 f 13 6s 2 2 F 7 / 2 metastable level 
of Yb + to be very sensitive to the Lorentz-violating effects. 
Indeed, we find that the LLI sensitive matrix element 
for this state is over an order of magnitude larger than for 
the 4f u 5d states. The results of the calculations are sum¬ 
marized in Table[l] where we also list T (2) matrix elements 
for nd states of Ca + , Ba + , and Yb + , which all have similar 
values. 


TABLE 1: Reduced matrix elements of T (2) operator in 
Ca + , Ba + , Yb + ions in atomic units. Ca + values are from 
Ref. Q. 


Ion 

State 

<7||r <2) ||7> 

Ca + 

3d 2 D}/2 

7.09(12) 


3d 2 D 5/ 2 

9.25(15) 

Ba + 

5 d 2 Dy2 

6.83 


5d 2 D 5/ 2 

8.65 

Yb + 

4f U 5d 2 Dy 2 

9.96 


4/ 14 5d 2 Z) 5/ 2 

12.08 


4/ 13 6i 2 2 F 712 

-135.2 


The calculations of the T (2) reduced matrix elements 
for the nd states were carried out using the all-order (lin¬ 
earized coupled-cluster) method l20l as described in the 
Methods section. 

In the case of Ca + ion, the frequency shift between the 
m = 5/2 and m = 1/2 3 d 2 D $/2 states is given by (7): 

\ (E mj= 5/2 - E mj= 1 / 2 ) = -4.45(9) x 10 15 Hz • C< 2) . (3) 

In the case of the 4 f 13 6s 2 2 / 7 7/2 state in Yb + , the frequency 
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FIG. 3: Relevant energy levels of Yb + with the hyperfine 
splitting shown for the 4 f 13 6s 2 2 / r 7/2 state of the 
171 Yb + isotope. 


shift between the m = 7/2 and m = 1/2 states is 

\ (E mj = 7/2 - £^= 1 / 2 ) = 6.14 x 10 16 Hz • C< 2) , (4) 

giving factor of 15 enhancement in comparison to Ca + due 
to the increased T (2) matrix element. A further significant 
advantage as compared to Ca + arises from much longer 
lifetime of this state as will be discussed below. 

We discuss now to detect potential LLI violations by 
measuring the energy difference between the mj = |7/2| 
and mj = |1/2| 4f 13 6s 2 2 F 7 /2 states in Yb + (see Fig. [ 2 ]). 

Typically, the main source of noise in such measure¬ 
ments is due to magnetic field fluctuations. In order to 
remove this noise, the work of Ref. Q used a superposi¬ 
tion of two ions with a state prepared in a decoherence-free 
subspace (DFS) |E I 22l . To implement the DFS technique 
with two trapped Yb + ions, one will monitor the phase 
evolution difference of the state IT 1 ) = -^=(| 1 / 2 , - 1 / 2 ) + 
|7/ 2 , -7/ 2 )), where \m \, m 2 ) represents the state with mj = 
mi and m 2 for the first and second ion, respectively, in 
the F 7/2 manifold. The target state IT 1 ) can be prepared 
by creating a product state, | X T P ^ = \ (|— 1 / 2 ) + |—7/2)) ( 8 ) 
(1+1/2) + |+7/2)), which dephases into a mixed state that 
contains IT 1 ) with 50% probability. A direct preparation of 
an entangled state would result in an increase of the signal- 
to-noise ratio by a factor of two, however, adds significant 
experimental complications. 

Starting from the ground state of Yb + , the state IT 1 ) = 

(| 1 /2, -1 /2)+17/2, -7/2)) can be prepared with 7r/2 and 
7r pulses using the coherent control developed for the im¬ 
plementation of the octupole 2 S 1/2 - 2 / 7 7/2 atomic clock 
mm Alternatively, one can also implement n/2 and 
7i pulses by driving Raman transitions via the 1 D[5/2]s/2 
state with 639 nm laser light (see Figure [3}. For the 
odd 171 Yb + isotope which has a nuclear spin of I - 1/2, 
one can prepare the target state IT 1 ) = \mp = 0,m/7 = 0) + 

| m F = 4, m F = -4) (F = 4) through the 2 F 7 / 2 (F = 3) state 
with radio-frequency pulses (see insert in Fig. [3}. 


The limits of the Lorentz-violating c^ v parameters ob¬ 
tained from measurements with Yb + are expected to im¬ 
prove greatly over the current limits due to two main fac¬ 
tors: sensitivity of the probe electronic state state and 
its lifetime. The experiment in Ref. (7) with 40 Ca + was 
performed with a Ramsey duration of 100 ms, with the 
main limitation coming from the heating of the ion mo¬ 
tion, which degrades the fidelity of the 7r/2 and n pulses 
during the state readout. This problem can be overcome 
with sympathetic cooling (see 120). However, the usable 
Ramsey duration is limited by the natural lifetime of the 
2 D 5/2 state approximately of 1.2 s to about 0.6 s. For the 
Ba + ion, the lifetime of the 2 D 5 / 2 state is approximately 
30 s (Ref. El), which allows for Ramsey durations of 
-15 s, potentially improving the measurement precision 
by V 15 /O .6 ~ 5 over 40 Ca + . 

Since the lifetime of the 2 F 7 /2 state of Yb + is in the 
range of several years |T8t the Ramsey duration of the pro¬ 
posed experiment with Yb + is not limited by spontaneous 
decay during the measurement. Instead, the coherence of 
the target state in the decoherence-free subspace is likely 
to be limited by the stability of the magnetic field gradi¬ 
ent in the vicinity of the ions. It has been shown that in a 
magnetically unshielded environment, ions can retain en¬ 
tanglement in the decoherence-free subspace for up to 30 s 
|[25l . Hence, in a magnetically shielded trap, we can expect 
to suppress decoherence due to magnetic field noise such 
that Ramsey durations for much longer times are realistic. 

The sensitivity of the frequency variation can be esti¬ 
mated from the quantum projection noise of the measure¬ 
ment: 

°7 = — 1 —~ ■ 2- (5) 

4 nAyfT yfr 

where A is the amplitude of the Ramsey oscillation signal, 
T is the Ramsey duration and r is the measurement time. 
Assuming that the target state IT 1 ) = -^(| 1/2, -1/2) + 

|7/2, -7/2)) for 2 F 7 / 2 is prepared through the creation of 
a mixed state as demonstrated in Refs. go, then the 
amplitude of the Ramsey oscillation is A = 0.5. With 
the Ramsey duration of 60 s, the estimated sensitivity is 
- 0.02/ V'r/Hz. Together with the improved sensitivity of 
the atomic state to Lorentz-violation as given in Eq. 0 . 
we estimate that the parameter C^ 2) will be bounded at the 
level of 1.7 x 10 -19 / V^/sec for two ions, which is an im¬ 
provement of - 2,400 times over the measurement with 
40 Ca + (Ref. (7)). Use of the pure entangled two-ion states 
would improve these bounds by an additional factor of 2. 

At this level of sensitivity, we expect that the c jk param¬ 
eters will be bounded at the level of 3x 10 -22 and 1.5x 10 -23 
for a day-long and year-long measurement, respectively, 
for two entangled ions (A = 1). The cjj terms, which 
are sensitive to the velocity of the laboratory frame, are 
suppressed by a factor of 10 -4 due to the relatively small 
velocity of the Earth with respect to the cosmic microwave 
background as compared to the speed of light. With a year¬ 
long measurement, one expects to be sensitive to the c F j 
parameters at the level of 1.5 x 10 -19 , a factor of 6,000 im¬ 
provement over the current best limits (6). In the electron- 
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photon sector, this would for the first time, probe all ten¬ 
sor Lorentz-violating elements below the ratio between the 
electroweak and Planck energy scales (2). 

Systematics which can directly affect the energy of our 
target states are magnetic and electric field fluctuations in 
the vicinity of the ions. Both fields can be measured di¬ 
rectly using the ions themselves as a probe. One can then 
use the knowledge of these fields to correct for the ef¬ 
fects of the differential linear Zeeman shift between the 
two ions due to a finite magnetic field gradient across the 
ion string, the quadratic Zeeman shift and the Stark shift 
from the oscillating electric field from the trap. A change 
in the environmental temperature during the measurement 
can also affect the energy shift through the black-body ra¬ 
diation (BBR). The BBR shift between the two Zeeman 
states (mj = 1/2 and nij = 7/2) of the F 7 / 2 manifold de¬ 
pends on the anisotropy of the BBR source. By assuming 
that the spatial temperature variation of the environment is 
less than 0.1 K, we estimate that in the worst case scenario 
of a maximal anisotropy of 100%, the differential BBR 
shift between the two Zeeman states is ~ 1 //Hz which 
affects the sensitivity of cjk at the 6 x 10 -24 level. 

Our calculations also identified general rules for the en¬ 
hancement of the reduced matrix elements of the T (2) op¬ 
erator. We find that the only parameter that significantly 
affects the matrix elements leading to the tensor LV sensi¬ 
tivity is the deeper localization of the probe electron. We 
find (if/\r\if/) of ~0.8 a.u or below for the corresponding 
electron to be a good indicator of the large value of T (2) 
matrix element. This condition is satisfied for the 4/ hole 
states, such as Yb + state considered here, or for highly- 
charged ions with nf valence electrons and degree of ion¬ 
ization - 15 [26). Such systems with configurations con¬ 
taining two 4/ electrons of two 4/ holes can have extra 
factor of two enhancement, and higher degree of ioniza¬ 
tion leads to further increase in the T (2) matrix element 
values. Other lanthanide ions, for example the Tm + ion 
with 4 f 13 6s ground state, may provide other potential can¬ 
didates for LV tests. We also considered the Th 3+ ion due 
to its 5 / 5/2 ground states and demonstrated laser cooling 
122, but we find that the corresponding T (2) matrix ele¬ 
ment is more that factor of 3 smaller than the one in Yb + . 
The 5/ orbital here is a valence orbital rather than a hole 
in the filled shell; moreover the 4/ electron is generally 
deeper localized than the 5/ one. 

In summary, we identify Yb + ions as an excellent sys¬ 
tem to test Lorentz symmetry in the electron-photon sec¬ 
tor. We also develop general rules for LV sensitivities in 
atomic systems. Systematic effects due to local magnetic 
field and electric field can be corrected for by using the 
ions themselves as an independent probes. The ultimate 
sensitivity of such experiments is likely to be limited by 
the differential BBR shift between the two Zeeman states, 
which we estimate to be approximately 10 times less than 
the projected sensitivity. 

Due to increased sensitivity to LLI-violation and life¬ 
time of the metastable F 7/2 state of Yb + , we estimate that 
experiments can reach the sensitivities of 1.5 x 10 -23 for 
the cjk parameters, more than 10 5 times stringent than the 


current best limits (7). Moreover, the projected sensitivity 
of the c T j parameters will be at the level of 1.5 x 10 -19 , be¬ 
low the ratio between the electroweak and Planck energy 
scales. 

It has been conjectured that Lorentz symmetry may be 
violated in string theories with the LV effects being sup¬ 
pressed by some power of m ew IM p i = 2 x 10 -17 , the ratio 
between the electroweak scale and the natural (Planck) 
energy scale for strings 0. If the Lorentz violation is 
not observed in the proposed here Yb + experiment, it 
will show that Lorentz violation in the photon-electron 
sector does not arise at this first-order level in strings, or 
a cancellation is present between photon and electron LV 
effects to nullify the combined result. With the proton and 
neutron tests already over the 0(1) bound (SHE), such 
result would demonstrate that if physics at the Planck scale 
violates Lorentz invariance, it is more than minimally 
suppressed at low energies for normal matter. 
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METHODS 


The calculations of the energy shift due to LLI violation 
reduces to the calculation of the expectation value of the 
Hamiltonian in Eq. The matrix element of (/m/|p 2 - 
3p\\Jmj) is expressed through the reduced matrix element 
of the r (2) operator using the Wigner-Eckart theorem 


[ J 2 3 ]</||r (2) ||7>. 

V -mj 0 mj ) 


Using the algebraic expression for the 3y-symbol, we ar¬ 
rive at the following expression for the matrix element of 
(Jmj |p 2 - 3p\\Jmj ) operator: 


(Jmj\T^\Jmj) = 


-J(J + 1) + 3m 2 , 


V(2/ + 3)(7+ l)(2J + l)J(2J-l) 
x </||r (2) ||/>. (6) 


The all-order method gave very accurate values of the 
3 dj lifetimes EH and quadrupole moments lf29ll in Ca + . 
In the all-order method, single, double, and partial triple 
excitations of Dirac-Hartree-Fock wave functions are in¬ 
cluded to all orders of perturbation theory. We refer the 
reader to review Ref. (20l for the description of the all¬ 
order method and its applications. The all-order results are 
accurate to about 1.5% a The calculations were carried 
out with both nonrelativistic and relativistic operators; the 
differences were found to be negligible at the present level 
of accuracy. The results are listed in Table [2] 

The all-order method was designed to work for mono¬ 
valent systems and is not applicable to the calculation of 
4/ 13 6 s 2 2 F 7/2 properties due to the electronic configuration 
of this state that has a hole in the 4/ shell. We used both a 
one-electron Dirac-Fock calculation and a large-scale con¬ 
figuration interaction calculation for this state and find ex¬ 
cellent agreement of both approaches. The 15-electron 
configuration interaction calculation follows the approach 
described in Ref. 130]. Briefly, we start from the solution 
of the Dirac-Fock equations and carry out the initial self- 
consistency procedure for the [ls 2 ,...,4/ 14 , 6 p] configura¬ 
tion. The 6s orbital was constructed for the 4f l3 6s 2 con¬ 
figuration, and the 5d^ 12 , 5/2 orbitals were constructed for 
the 4f l3 5d6s configuration. The basis set used in the Cl 
calculations included virtual orbitals up to 8 s, 8 p, Id , 7/, 
and 5 g. The virtual orbitals were constructed as described 
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TABLE 2: Reduced matrix elements of T (2) operator in 
Ca + , Ba + , Yb + ions in atomic units. Ca + values are from 
Ref. El • Method of calculations is listed in the column 
“Method”. LCCSD is linearised coupled-cluster method 
with single and double excitations, DF is Dirac-Fock 
method, RPA is random-phase approximation, and Cl is 
configuration iteration. We list the included 
configurations for the 15-electron Cl calculation. 


Ion 

State 

Method 

(J\\T (2) \\J) 

Ca + 

m 2 d 3/2 

LCCSD 

7.09(12) 


3 d 2 D$/2 

LCCSD 

9.25(15) 

Ba + 

5d 2 D 3/2 

LCCSD 

6.83 


5 d 2 D$/2 

LCCSD 

8.65 

Yb + 

4f H 5d 2 D 3/2 

LCCSD 

9.96 


4f H 5d 2 £>5/2 

DF 

7.23 



Cl Isp6df5g 

11.6 



LCCSD 

12.08 


4f l3 6s 2 2 F 7/2 

DF 

-145 



Cl 6sp5df 

-143.8 



CI+RPA 6sp5df 

-135.1 



CI+RPA Isp6df5g 

-135.6 



CI+RPA 8sp6df5g 

-135.2 


in Refs. [31, 32). The configuration space was formed by 
allowing single and double excitations for the odd-parity 
states from the 4 f u 6p, 4 f 13 6 s 2 and 4/ 13 5d6s configura¬ 
tions to virtual orbitals of the basis set listed above. 


We have verified the convergence of the Cl by carrying 
out three calculations with an increasing set of configura¬ 
tions functions: (1) including single and double (SD) ex¬ 
citations to the 6 s, 6 p, 5d, and 5/ orbitals (we designate 
it [ 6sp5df ]), (2) adding excitations to the 7s, Ip , and 6 d , 
6 / and 5 g orbitals [7 sp6df 5 g], and (3) also adding 8 s, 8 p, 
Id , and 7/ orbitals [8spldf5g]. The last configuration set 
results in rather lengthy calculations with ~ 2 300 000 de¬ 
terminants. 


As demonstrated in Table[2j the number of included con¬ 
figurations has only a negligible effect on the T (2) operator 
of the 4/ 13 6s 2 2 F 7/2 state. The Cl number is in agreement 
with simple DF calculation. Inclusion of the random-phase 
approximation changes this value by 6 %. 







